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On the light-bending model of X-ray variability of MCG— 6-30-15 
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Abstract 

Wc apply the light bending model of X-ray variability to Suzaku data of the Seyfert 1 galaxy MCG-6- 
30-15. We analyze the energy dependence of the root mean square (rms) variability, and discuss conditions 
necessary for the model to explain the characteristic decrease of the source variability around 5-8 keV. A 
model, where the X-ray source moves radially rather than vertically close to the disk surface, can indeed 
reproduce the reduced variability near the energy of the Fe line, although the formal fit quality is poor. 
The model then predicts the energy spectra, which can be compared to observational data. The spectra 
are strongly reflection dominated, and do not provide a good fit to Suzaku spectral data of the source. 
The inconsistency of this result with some previous claims can be traced to our using data in a broader 
energy band, where effects of warm absorber in the spectrum cannot be neglected. 
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1. Introduction 

The Seyfert 1 galaxy MCG~6-30-15 shows the best ex- 
ample of a very broad Fe Kq, fluorescent line in its X-ray 
spectrum (Tanaka et al. 1995; see Reynolds & Nowak 2003 
for review). According to numerous studies (e.g., Miniutti 
et al. 2007 and references therein) the line is produced 
well inside the marginally stable orbit for a Schwarzchild 
black hole, requiring a (maximally) rotating black hole. 
Uncertainties remain, though, as to the robustness of the 
broadness of the line, since the source has a very complex 
ionized absorber (e.g., Otani et al. 1996; Lee et al. 2001) 
and the details of the line modelling (e.g. its width) are 
rather sensitive to the description of the absorber (Miller, 
Turner & Reeves 2008; Miyakawa, Ebisawa & Inoue 2010). 

At the same time, the variability properties of the line 
are rather incompatible with a simple geometrical ideas of 
the central accretion flow. The line, and the Compton 
reflection component associated with it appear to be much 
less variable than the continuum producing the repro- 
cessed component. This was pointed out by Inoue & 
Matsumoto (2001, 2003), and later confirmed by many 
authors (e.g., Fabian et al. 2002). The line does not follow 
the continuum variability in a simple way, although some 
changes in the line profile/flux are observed (Iwasawa et 
al. 1996). 

In simple geometrical scenarios the line and reflected 
continuum should respond to rapid variations of the pri- 
mary continuum driving the reprocessed components. An 
explicit decomposition of the MCG-6-30-15 spectra onto 
components which can be interpreted as variable and con- 
stant one was done by Miller et al. (2008) using the 



Principal Component Analysis method. The constant 
component appears to have a shape of a mildly ionized 
reflection, while the variable one corresponds to a broad- 
band continuum with ionized absorption. 

A way of reconciling these two properties of the Kq 
line was the so called "light bending model" (LBM) by 
Miniutti & Fabian (2004). The idea is that, for a source 
close to the central black hole, light bending and other 
relativistic effects will affect much more the observed pri- 
mary emission than the reflected component. Thus, as 
the position of the source changes, the observed flux of 
the primary source will vary more than the flux of the re- 
flected component. The model gained popularity despite 
it being rather ad hoc and having no obvious relation to 
the whole body of knowledge on X-ray variability of ac- 
creting black hole systems (see M'^Hardy 2010 for a recent 
review and references). 

The properties of the LBM were reanalyzed by 
Niedzwiecki & Zycki (2008, hereafter NZ08), where it was 
concluded that in the original formulation the model ac- 
tually failed to explain the data (sec also Niedzwiecki & 
Miyakawa 2010, hereafter NMIO). NZ08 suggested an al- 
ternative geometrical scenario, where the source moves 
radially (rather than vertically, as in the original formu- 
lation), very close to the disc surface. Then, light bend- 
ing to the disc plane, an effect specific to the (extreme) 
Kerr metric causes the blue peak of the line to be approxi- 
mately constant. The best variant of the model to explain 
the lack of variability of the line was that assuming the in- 
trinsic source luminosity followed the Keplerian accretion 
disc emissivity (see fig. 7 in NZ08). Here, for radii < 3i?g, 
where Rg = GM/c^, the line flux is almost constant, while 
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the primary source flux changes by a factor of « 7. 

In the following paper NMIO extended the model calcu- 
lations to include the reflected continuum. The full model 
was then applied the Suzaku rms(i?) spectra of MCG-6- 
30-15. The rms(£') spectrum is defined as follows, 



rms(-B) 



1 



{F{E))' 



N-1 



where the sum is over N contributions to flux F{E) at en- 
ergy i?, while {F{E)) is the mean value of the flux. The 
rms(i5) is thus a measure of the source variability, as a 
function of energy, E. The rms(i?) spectra could indeed 
be modelled assuming that the source is located very close 
to the black hole, R = 1.6-3i?g. The rms(£^) spectra com- 
puted assuming that spectra from 3-5 random radii within 
that range contribute to the emission could indeed match 
the shape of the observed rms(£'). However, to match the 
amplitude of the rms(£') the intrinsic source normaliza- 
tion had to be assumed to vary. The intrinsic variability 
deprives the model its major attracting feature, namely 
the explanation of the constancy of the reflected compo- 
nent. On the other hand it may go some way towards 
explaining other key features of observed X-ray variabil- 
ity of Active Galactic Nuclei: the power spectrum, hard 
X-ray time lags, etc. 

In this paper we investigate the model properties in 
some more detail. Our main purpose is to compute the 
energy spectra from the model and compare them with 
the Suzaku spectral data of MCG-6-30-15. 

2. The light bending model 

Our computations of the model are described in detail 
in NZ08 and NMIO, following the original formulation by 
Miniutti & Fabian (2004) . Here we give only a brief sum- 
mary of its assumptions and properties. We assume a 
maximally rotating Kerr black hole and we consider the 
geometry, where the X-ray source moves radially, very 
close to the disk surface. The motion is along the line of 
constant polar angle, 9, with H = 0.07R, where H is the 
height of the source. The intrinsic luminosity of the source 
follows the Keplerian accretion disk emissivity, L^'^{R), 
(Page & Thorne 1974; model S^^ in NZG8), where 



3GMM 
87ri?3 



where r = R/ Rg, and i?fl(r,a) is the relativistic correction 
factor, containing the zero-inner torque condition and de- 
pendent on radial distance r and black hole angular mo- 
mentum a [see Krolik 1999 for explicit form of i?R(r, a)] 

Then NZ08 show (see fig. 7 there) that, for R<3Rg the 
observed fiux of the Fc line is practically constant while 
the observed flux of the primary continuum varies by a 
factor of « 7. We illustrate this result in Fig. 1, which 
shows the model spectra for different radii < 3i?g. The 
flux of the reflected component (uniquely related to the 
flux of the Fe Kq, line for a given primary spectral shape 
and reflector ionization) is practically constant, while the 
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Fig. 1. The light bending model energy spectra from 
the model of NZ08, i.e. assuming that the intrin- 

sic source luminosity follows the accretion disk dissipation 
rate (for a=0.998) and the source can move along the line 
of constant polar coordinate with its height H = 0.07-R. 
Observation (inclination) angle is i = 35° and the source 
radial positions (corresponding to different curves) are 
1.6,1.8,2,2.2,2.4,2.6,2.8,3iig. The observed level of the direct 
source emission (represented by the approximately horizontal 
lines) depends then strongly on the source position, while the 
reflected component (and, in consequence the Fe line) is al- 
most constant (see also fig. 7c in NZ08). 



primary emission changes significantly. In the following 
computations we use the model in the version extended 
by NMIO, where the radial location of the X-ray source 
is generated from a distribution P{r) oc , with 6 = —1.5 
(sources concentrated towards the black hole) , and the in- 
trinsic luminosity is randomized by being generated from a 
gaussian distribution centered at L^"^(i?), with standard 
deviation cr(i?) = 0.2L^"^(i?) (note that the randomiza- 
tion of luminosity is not included in the spectra plotted 
in Fig. 1). 

We note here that the model predicts spectra which are 
strongly reflection dominated. This is a necessary con- 
sequence of the basic assumptions of the light bending 
model. The major problem then the model will face is 
whether such reflection dominated spectra can be com- 
patible with MCG-6-30-15 data. We note that Fabian et 
al. (2002) did suggest that the X-ray spectra of Seyfert 1 
galaxies may be reflection dominated but no fits to data 
of MCG-6-30-15 were presented. Since then, such reflec- 
tion dominated spectra have been claimed to be good fits 
to data from a number of objects, including MCG-6-30- 
15. Most recent analysis of this source was published by 
Miniutti et al. (2007), where Suzaku data were analyzed. 
The large value of the amplitude of reflection inferred from 
the fits, il/ {2tt) w4, is claimed to be consistent with LBM. 
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Fig. 2. The rms(_E) dependence from Suzaku data of 
MCG— 6-30-15, together with model curves. Data points show 
rms(_E) for time bin of 16384 sec (upper set) and 131072 sec 
(lower set). Each pair of curves shows a ±1(7 range of rms{E) 
from an ensemble of realizations of the model. Numeric labels 
show riR - the number of radii contributing to each spectrum 
(see text, Sec 3, for details). 

3. Modelling the energy dependence of the r.m.s 

Results of NZ08 imply that a sequence of spectra from 
R<3Rg will produce a rms(ii') with a dip around the line 
energy, since the amplitude of reflected continuum is the 
same for all radii (Fig. 1). 

Physically, we can expect that a spectrum from a given 
time interval of observation, T, will consist of contribu- 
tions from a number of radii, the number depending on 
T. We treat the number of radii, riR, as a parameter of 
the model, since the exact physical processes responsible 
for generating active regions are not known. The rms(£') 
spectrum is thus computed generating values of radius 
from the 1.6 Rg < R < 3i?g range, adding energy spectra 
from these radii to form one basic spectrum, and then re- 
peating the process N times, which corresponds simply 
to dividing the whole observation into N intervals. The 
rms(£') is computed from these N spectra. 

A single rms(£') spectrum is just one realization of a 
random process. Repeating the above procedure would 
produce different rms(i!^). To estimate the range of possi- 
ble ims{E) behavior we repeat the procedure 1000 times, 
and plot the resulting ±la range of ims{E) in Fig. 2. 
Since we do not perform formal fitting [x^ minimization) 
of the model to the ims{E) data, the plot helps to es- 
timate the optimum value of a parameter, in this case 
riR. Specifically, rms(£') computed with rip ~ 2-3 gives 
the best match to the shorter time scale (16 ksec) ob- 
served rms(£^), while tir = 4-5 might be appropriate for 
the longer time scale ims{E), although in the latter case 
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Fig. 3. Contributions to from fitting the spectral model 
predicted by the LBM matching the rms(i?) relation, to the 
time averaged Suzaku data of MCG-6-30-15. The fit is not 
satisfactory, with xi ~ The model spectrum is strongly 
reflection dominated and extremely relativistically smeared, 
since the source radial position R, < 3 Rg . 

the slope of the model ims{E) is rather flatter than that 
of the data. The normalization of the model rms(£^) spec- 
trum is not adjusted to best match the data. This is dif- 
ferent than the usual practice with energy spectral fitting, 
where the model amplitude is adjusted to give lowest pos- 
sible ■ Therefore the formal values computed from 
these rms(£') models are rather high, the lowest values of 
per data point being « 2. 

4. Modelling the energy spectra 

4-.1. Energy spectra predicted by the light bending model 

We now present the energy spectra implied by the fits 
of the LBM to the ims{E) data and attempt to compare 
these with the spectral data of MCG-6-30-15. 

Each realization of the rms(-E) yields certain mixture 
of energy spectra from a number of radii, with relative 
weights resulting from the normalizations generated in the 
process (depending on the radius and the random ampli- 
tude around the L^'^(r) value). The total spectrum is 
then a sum of all these spectra. The spectra are strongly 
reflection dominated (Fig 1) with the amplitude of reflec- 
tion, r2/(27r) reaching a factor of « 5. 

We fit the January 2006 Suzaku XIS and PIN time av- 
eraged data in the range 1-40 keV. Importantly, in this 
energy range the effects of warm absorber cannot be ne- 
glected, so we fit a model which includes description of 
the absorber. Specifically, we apply a warm absorber 
model constructed from XSTAR calculations (Kallman 
et al. 2004). We find that a double model (low and high 
ionization) is necessary and sufficient to model the data 
(Sec 4.2; see also Miyakawa et al. 2009 for details of data 
reduction and warm absorber modelling). 
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We do preliminary fitting of the spectra from the NZ08 
LBM model to the data. The source's radial position is 
the most important fitting parameter. Computations of 
the reflected component were done only for cold plasma 

= 0) and for Solar abundances of metals (constructing a 
grid of models is not practical due to computing resources 
required), therefore the model is not as general as other 
models available in XSPEC. The fits are generally poor, 
with reduced = 1-8 — 1.9. Strongest residuals appear 
below 5 keV. 

We construct a more general reflection model by us- 
ing the results of computations by Ross & Fabian (2005). 
Their model (table model ref lion in XSPEC) includes 
all relevant atomic physics, including the Fe line emission 
(although it is not valid for low ionization, ^ < 30). It 
was also used in many previous papers (e.g., Miniutti et 
al. 2007). We now apply relativistic effects by convolv- 
ing the ref lion spectrum with the transfer function from 
our LBM calculations for a given radius. Then, we add 
a power law continuum normalized so that the amplitude 
of the reflected component corresponds to that from the 
LBM computations. To achieve this, we use the ratio of 
the numbers of photons in the primary and reflected con- 
tinua above 12 kcV (the limit is important to exclude the 
effects of ionization) . Finally, we sum the spectra over the 
set of radii. The resulting model is thus parameterized by 
primary power law slope, F, ionization parameter, ^, Fe 
abundance and the inclination angle. The amplitude of 
reflection and relativistic smearing effects (i.e. the radial 
emissivity) result form the properties of the LBM and are 
thus fixed for a given set of radii. 

This spectral model gives somewhat better fit, but it 
is still not satisfactory, with xt ~ 1-5 (Fig 3). Strong en- 
hancement of reflection drives the primary continuum to 
be steep, F w 2.4, which seems incompatible with the slope 
required for the lower energy range (Fig 4) and, in conse- 
quence, significant deviations in are present at _E < 6 
keV. The best fits are obtained for inclination of 40° and 
Fe abundance equal to Solar. The quality of spectral fits 
is basically independent of the model realization, contrary 
to the quality of the rms(i?) fits. 

4-2. Time averaged spectrum of MCG-6-30-15 

According to some earlier studies, the reflection dom- 
inated spectra provide good description of the MCG-6- 
30-15 data, which fact has been used as a supporting 
argument for the applicability of the LBM. In particu- 
lar, Miniutti et al. (2007; see also references therein) fit- 
ted the January 2006 Suzaku data with a model com- 
prising cutoff power law continuum and ionized reflection 
(XSPEC ref lion model of Ross & Fabian (2005) with 
relativistic smearing modelled with extreme Kerr metric 
transfer function of Laor (1991; XSPEC kdblur2 model). 
Using this model Miniutti et al. obtain their best fit with 
strongly enhanced reflection, f2/(27r) « 4. The reflected 
component is strongly smeared, with inner disk radius of 
1.7 -Rg and illumination emissivity following very steep, 
r~'*'^, dependence up to the break radius of 6Rg, beyond 
which the emissivity is flatter, r~^'^. 



Confidence contours: Chi-Squared 
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Fig. 4. The contours in the f2/(27r)-r plane for model fits 
to Suzaku data of MCG-6-30-15 in two cases: data above 3 
keV, no warm absorber model (big contours) and data above 
1 keV with a warm absorber model included (small contours). 
There is a clear and significant dependence of the amplitude 
of reflection on the data range/model components used. 

In order to understand then why our results disagree 
with those earlier ones, we compare the results of fitting 
the same data with the same reflection model (i.e., the rel- 
ativistically smeared ref lion model) in two cases: (1) us- 
ing only data above 3 keV, assuming that warm absorber 
does not affect the spectrum (as in Miniutti et al. 2007), 
and (2) using the data above 1 kcV and attempting to 
describe the warm absorber. The warm absorber model, 
computed with XSTAR is prepared as a table model on 
a grid of 20x20 values of column density, N^, and ion- 
ization parameter, ^ (Miyakawa et al. 2009). We use two 
tables (lower and higher ionization), having checked that 
a single-^ model does not provide a good fit, while adding 
a third table docs not improve the fit. We assume that 
the radial emissivity for line profile computation is the 
same as in Miniutti et al. (2007) and we keep it fixed, 
but we let the inner disk radius and inclination to be free. 
We keep the outer disk radius fixed at 400 i?g. The Fe 
abundance is also free to fit, as is obviously the overall 
reflection amplitude. 

We reproduce the results of Miniutti et al. (2007) in the 
first case, where only data above 3 keV are used. However, 
with broader energy band data the best-fit results are very 
different. The spectral index is smaller (harder spectrum) 
and the reflection amplitude is much smaller, in fact less 
than 1. In Fig. 4 we show the x^ contours in the spectral 
index - reflection amplitude plane, for both considered 
cases. The inner disk radius, characterizing the broaden- 
ing of the line is comparable in both cases. 

Our results mean that, to state it conservatively, the de- 
rived parameters of the reflection component are strongly 
model dependent. Fits in broader energy band yield 
harder spectrum and, in consequence, much smaller am- 
plitude of reflection than fits performed in the E > 3keV 
energy band. This is a problem for LBM since it neces- 
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sarily predicts reflection dominated spectra. 

It is not our goal here to construct more complex models 
and to study the dependence of the reflection parameters 
of the way how absorption is modelled. It seems that the 
solutions of Miller et al. (2008) and Miniutti et al. (2007) 
cover the whole broad range of possibilities. We simply 
want to emphasize that, in the context of our considera- 
tions, the support for the LBM depends on modelling the 
complex spectrum of MCG-6-30-15. 

5. Discussion and conclusions 

The light-bending model is an attempt to reconcile the 
extreme width of the Fe line with its lack of variability, 
first reported in MCG-6-30-15. It invokes extreme rela- 
tivistic effects to explain weaker variability of the repro- 
cessed spectral components compared to the primary com- 
ponents. Following the original formulation of the model 
by Miniutti & Fabian (2004) we reconsidered the model 
in greater details (NZ08, NMIO) pointing out some inac- 
curacies in the original version. Reduced variability of the 
Fe line compared to primary continuum is indeed possible, 
but for a radially rather than vertically moving source. 

Quantitatively, the model is able to reproduce the 
rms(£') relation, in particular the drop near the Fe line 
energy, but only if some intrinsic variability of the X-ray 
source is assumed (NMIO). That is not necessarily a bad 
thing since it might go some way towards reproducing 
the usual characteristics of X-ray variability (e.g., power 
spectra), which the original model does not address. One 
possible test of the model is then to construct the energy 
spectra and compare these to the data. The model spectra 
are obviously strongly reflection-dominated and relativis- 
tically smeared. 

Wc find that the refiection-dominated spectra do not 
fit the data, contrary to previous studies (Miniutti et al. 
2007 and references therein). There is a significant de- 
pendence of the best-fit models on the energy range of 
the fitted data. The effects of warm absorber in the spec- 
trum are not limited to below 3 keV, and so using the data 
below this energy and including a warm absorber model 
gives very different results compared to fits without the 
warm absorber and data above 3 keV. In this respect our 
results can be located between the results of Miller et al. 
(2008) and Miniutti et al. (2007), which can be thought of 
as absorption-dominated and reflection-dominated solu- 
tions, respectively. We note that the dependence of reflec- 
tion parameters on details of the warm absorber descrip- 
tion applies also to the width of the Fe Kq line. While 
this is usually claimed to be extremely broad, data mod- 
elling including complex warm absorber models (Miller et 
al. 2008; Miyakawa et al. 2010) yields good descriptions of 
the broad band data of MCG-6-30-15 without extremely 
broad Fe line. 

MCG-6-30-15 is not the only source showing certain 
degeneracy between reflection and absorption in their X- 
ray spectra. Similar situation occurs in the X-ray bi- 
nary XTE J1650-500 where Miniutti, Fabian & Miller 
(2004) presented fits of models with extremely broad Fe 



line, while Gicrlihski & Done (2006) demonstrated that 
the Fe line docs not need to be so broad, if complex 
absorption and rcfiection models are applied. Another 
spectral feature where this kind of ambiguity is observed 
are the so called soft X-ray excesses observed in many 
Active Galactic Nuclei. These can be modelled equally 
well by absorption and reflection effects (Sobolewska & 
Done 2007). The situation may appear somewhat frus- 
trating, appearing despite years of data collecting and ef- 
forts to understand them. More importantly though, the 
implications of the two models are very different and far- 
reaching, since they touch the base of our interest in X-ray 
astronomy. Studying effects of extreme gravity near com- 
pact objects was and still is one the main driving forces 
for developing technologies for X-ray observations. 

Summarizing, we fit the spectra predicted by the LBM 
of NZ08 to Suzaku data of MCG-6-30-15 and find that 
the model does not fit the data, if the fits use data in the 
energy range where warm absorber effects are important. 
Detailed results depend on models of warm absorber, with 
more complex warm absorber models implying less ex- 
treme reflection. 
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